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The  cumulative  distribution  of  the  number  of  secondary  electrons  in  a 
single-stage  photomultiplier  is  calculated  by  numerically  integrating  the  in¬ 
version  Integral  for  its  probability  generating  function  along  a  suitably 
chosen  contour.  A  residue  series  applicable  in  certain  cases  is  also  pre¬ 
sented.  Saddlepoint  approximations  to  the  contour  Integral  are  described, 
vhlch  are  the  more  accurate,  the  greater  the  numbers  of  secondaries.  Recur¬ 
rent  relations  are  developed  for  computing  values  of  the  distribution  for 
purposes  of  comparison.  Computation  of  the  Neyman  Type-A  distribution  is 
treated  as  a  limiting  case. 
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I.  Particle  Multiplication 


(a)  Introduction 

The  particle-counting  distributions  with  whose  computation  this  paper  is 
concerned  are  exemplified  by  that  of  the  output  of  a  photomultiplier  with  a 
single  stage  of  multiplication.  In  this  device  primary  photoelectrons,  driven 
out  by  incident  light,  are  accelerated  in  an  electric  field  and  impinge  on  a 
surface  from  which  they  eject  secondary  electrons.  Let  be  the  number  of 
secondary  electrons  ejected  by  the  j-th  primary  electron.  Then  if  k  primary 
electrons  strike  the  surface  during  a  fixed  Interval  (0,  T),  the  total  number 
of  secondary  electrons  is 

k 

n-J^Xj  (1.1) 

J-1 

The  nuiid>er  k  of  primary  electrons  is  a  random  variable  with  probability  gen¬ 
erating  function  (p.g.f.) 

f(z)  (1.2) 

k-0 

where  is  the  probability  that  k  primary  photoelectrons  strike  the  surface. 

Let  the  nuiid)ers  x^  of  secondary  electrons  be  Independent  and  identically  dls- 

(s) 

trlbuted  random  variables  with  probabilities  Pr  (x^  *  m)  -  p^  and  p.g.f. 

«o 

g(*)  *  •  (1.3) 

«i“0 

Then'  the  p.g.f.  of  the  total  number  n  of  secondary  electrons  counted  during 
<0,  T)  is 

m 

h(z) 
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where  is  the  probability  that  the  total  number  of  secondaries  equals  n  [1]. 

Although  the  probability  distribution  {p  }  may  be  of  Interest  In  comparing 

n 

theory  with  experiment,  the  cumulative  probabilities 

n-1 


m«0 


1,  2,  ... 


Qn  -  0. 


and  their  complements 


(1.5) 


S' 

i 


i 


■  * 


(1.6) 


m»n 


ere  generally  more  useful  because  they  characterize  the  perfoinnance  of  devices 
In  which  some  action  Is  triggered  when  the  number  n  of  particles  exceeds  a  cer¬ 
tain  bias  level.  False-alarm  and  detection  probabilities  In  optical-  and 
particle-detection  systems  and  error  probabilities  In  optical  communications 
are  directly  related  to  the  cumulative  distribution,  and  It  Is  the  computation 
of  this  distribution  that  will  be  studied  here.  We  call  and  the  "tall 
probabilities."  Their  generating  functions  are 


n®0 


H 


n*0 


(1.7) 


(1.8) 


In  terms  of  the  p.g.f.  h(z). 

Multiplicative  processes  such' as  this  occur  In  many  other  contexts  as  well. 
A  review  has  recently  been  given  by  Teich  (2],  who  mentions  applications  to 
biology,  medicine,  cosmology,  geophysics,  and  operations  research.  The  prob¬ 
abilities  p^  and  the  tall  probabilities  and  must  usually  be  determined 


-  I:-'.:- -I'v'-i  -I- 


by  numerical  methods,  but  the  conventional  ones  run  Into  difficulties  when  the 
nuBbers  n  and  their  expected  value  E(n)  are  large.  The  number  of  steps  needed 
to  compute  p^  generally  Increases  with  Increasing  n,  as  does  the  number  of  terms 
to  be  summed  In  evaltiating  .  Inaccuracy  In  the  Individual  values  of  p^  Intro¬ 
duces  large  relative  errors  Into  the  complementary  cumulative  probabilities 
a  1  -  qJ*  idien  n  much'  exceeds  E(n) . 

In  this  paper  we  shall  present  methods  for  computing  the  tall  probabilities 
that  become  the  more  efficient,  the  greater  the  numbers  n  and  E(n),  and  do  not 
require  computing  the  Individual  probabilities  p^  and  summing  them.  The  appli¬ 
cation  to  the  output  of  a  multiplicative  process  may  be  regarded  as  exempli¬ 
fying  methods  that  are  worthy  of  consideration  whenever  the  p.g.f  h(z)  of  some 
Integer-valued  random  variable  n  is  known  In  analytical  form  and  the  expected 

value  E(n)  and  the  numbers  n  for  which  Q  or  Is  wanted  are  large. 

n  n 

The  principal  method  Involves  suitably  deforming  the  contour  of  integration 
in  the  Inversion  integrals 

^  4k 

c“ 

Q+  -  f  _d2_  .  . 

'^n  /  2-1  27Ti»  U.IO) 

•'c 

idiere  C  and  C***  Initially  are  closed  contours  enclosing  the  origin,  but  no 
singularities  of  the  m.g.f.  h(z);  encloses  the  point  z  «  1,  C*  does  not; 
and  the  contours  are  traversed  counterclockwise.  The  Integration  along  the  de¬ 
formed  contour  Is  evaluated  by  the  trapezoidal  rule.  The  use  of  numerical  con¬ 
tour  Integration  to  evaluate  cumulative  distributions  of  continuous  random  vari¬ 
ables  was  treated  In  [3],  and  a  comment  on  that  paper  described  Its  application 
to  cumulative  distributions  of  Integer-valued  random  variables  [4]. 
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contour 


Besides  evaluating  the  tall  probabilities  and  numerical 

Integration,  we  shall  show  how  to  approximate  them  by  Isolating  the  contribu¬ 
tions  to  the  integrals  In  (1.9)  and  (1.10)  at  the  saddlepolnts  of  the  Integrand. 
Daniels  [5]  first  demonstrated  saddlepolnt  approximations  to  probability  dis¬ 
tributions  integer-valued  random  variables.  Related  to  this  method  is 

the  use  of  tilted  distributions,  utilized  by  Blackwell  and  Hodges  [6],  Bahadur 
and  Rao  [7],  Petrov  [8],  and  Bamdorff-Nlelsen  and  Cox  [9]  to  calculate  cumu¬ 
lative  distributions  of  sums  of  independent  random  variables;  see  also  Van  Trees 
[10].  A  different  saddle-point  approximation,  which  avoids  the  use  of  the 
error-function  integral,  was  utilized  for  tall  probabilities  of  integer-valued 
random  variables  in  [11].  In  the  present  problem  the  contributions  of  saddle- 
points  above  and  below  the  real  axis  in  the  z-plane  must  be  Included  in  addi¬ 
tion  to  that  of  the  principal  one  lying  on  the  real  axis.  We  shall  evaluate 
the  latter  by  means  of  a  uniform  asymptotic  expansion  [12] ;  the  contributions 
of  the  off-axis  saddlepolnts  are  evaluated  by  the  method  of  [11]. 

The  contour- Integral  method  and  its  approximations  apply  whenever  the  p.g.f. 
h(z)  of  the  output  distribution  is  known  in  analytical  form,  but  for  the  sake 
of  definiteness  we  shall  restrict  our  discussion  to  distributions  arising  when 
the  light  ejecting  the  primary  electrons  is  incoherent  light  with  various  simple 
spectral  densities,  and  the  number  of  secondary  electrons  ejected  by  each  prl- 
amry  electron  is  governed  by  a  Poisson  distribution.  A  by-product  of  our  study 
is  an  illustration  of  how  the  output  distribution  depends  on  the  product  of  the 
banihrldth  U  of  the  incident  light  and  the  duration  T  of  the  counting  Interval. 

We  first  review  methods  of  computing  the  probabilities  p^  that  follow  most 
directly  from  (1.4)  and  are  useful  when  the  numbers  n  and  E(n)  are  small,  in 
order  to  generate  accurate  values  of  those  probabilities  for  comparison  with 
the  contour- integration  and  other  methods  to  be  presented  later. 


(b)  General  Multiplicative  Processes 


From  the  p.g.f.  h(z)  the  probabilities  p^  can  be  determined  by 


Pq  • 


-L  d"h(z)|  -Id” 


n  n!  ^  n  I  ^  n!  ^  n 


f(g(2))  .  n  >  0. 

dz“  'z=0 


(1.11) 


dz“  'z=0 

By  means  of  a  formula  based  on  Bell  polynomials  [10,  p.l998],  the  probability 


p  can  be  written  as 
*^n 


'n-E 


k-1 


f,  H  ,  ,  nil, 
k  o,k*  * 


where 


dz  'z*g(0) 


g(0)  -  p, 


(s) 


and  the  H  ,  are  determined  by  the  recurrent  relation 
n,k  ^ 

®n+l,l  Pn+1* 


(1.12) 


(1.13) 


®n+l,k+l  “  5^  "  n+1^  Pn+l-m“m,k» 

m«K 


(I.IA) 


in  terms  of  the  distribution  of  the  number  of  secondaries  per  primary  photo¬ 
electron.  In  particular 

(1.15) 


n  .  «  (®>  /n! 


(c)  Poisson-Distributed  Secondaries 

When,  as  we  assume  henceforth,  the  distribution  of  the  number  of  secondaries 
per  primary  electron  has  the  Poisson  form  with  mean  G,  which  is  called  the  gain, 

pj*^-  g"  e’^/n!,  (1.16) 

the  probabilities  p  that  the  total  number  of  secondaries  equals  n  form  what 

II  « 


5  * 


is  called  a  compound  or  generalized  Poisson  distribution  [13]  and  are  given  by 

OD 

^n  "  2  V*^®/n! .  (1.17) 

k-0 

The  p.g.f.  of  the  output  distribution  is  now 

h(z)  -  f(e®^*"^^).  (1.18) 

It  follows  from  Problem  26,  p.  47  of  Riordan's  book  [14]  that 

-  g"  S(n.  k)  (1.19) 

In  terms  of  the  modified  Stirling  numbers  of  the  second  kind,  S(n,  k) ,  which 
obey  the  recurrent  relation 

S(l,  1)  -  1,  S(k.  n)  -  0.  k  >  n. 

S(n+1,  k)  -  [S(n.  k-l)+kS(n,  k)]/(n+l).  (1.20) 

In  terms  of  the  ordinary  Stirling  n\]inbers  of  the  second  kind  these  are  de¬ 
fined  by 

S(n,  k)  =  S(n,  k)/n!  (1.21) 

and  are  Introduced  in  order  to  avoid  overflow  in  machine  computation.  Recur¬ 
rent  relations  for  S(n,  k)  are  given  in  [14].  In  particular 

S(n,  1)  -  S(n,  n)  -  1/n!.  (1.22) 

As  a  result  the  probabilities  sought  are 


Pq  “ 


f^S(n.  k)e"''^  n  >  0, 


(1.23) 


k-1 


In  which  the  derivatives  f^^  of  the  p.g.f.  f^z)  of  the  primary  distribution  are 
-G 


This  enables  computation  of  the  probabilities  p^  by  a 


evaluated  at  z  ■  e 
strictly  finite  procedure  in  contrast  to  the  infinite  series  in  (1.17), 
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When  the  number  n  is  large,  however,  the  great  number  of  Iterated  com¬ 
putations  required  by  (1.21-.23)  may  introduce  substantial  errors  due  to 
rounding  off.  Greater  accuracy  can  be  achieved  for  large  n  by  utilizing  (1.17), 
replacing  n!  by  Stirling's  approximation  [15]  and  writing  it  as 

p^  -  [2Ti(n+l)]'^lB(n+l)r^ 

k-1 

where  ^ 

-  “p  (ik  -  ^ 

The  summation  was  halted  when  the  terms  in  (1.24)  became  insignificant. 


II.  Primaries  Ejected  by  Incoherent  Light 


(a)  Arbitrary  Spectral  Density 

The  incident  light  is  assumed  to  be  quaslmonochromatlc,  linearly  polarized. 
Incoherent  light  with  a  spectral  density  4(b));  the  angular  frequency  U)  is  mea¬ 
sured  from  the  central  angular  frequency  of  the  light.  Because  the  field  of 
the  light  is  a  Gaussian  random  process,  the  probability  generating  function 
of  the  distribution  of  primary  pfiotoelectrons  is 


w 

f(*)  “  J~[  tl-N  X  (^-1)1“^, 


where  is  the  mean  number  of  primary  electrons  [16]. 
values  of  the  Integral  equation 


(2.1) 


The  are  the  eigen- 


/■ 


Xp(t)  -  I  ^(t-8)p(s)  ds, 

^0 


(2.2) 


whose  kernel 


♦(T) 


I  ♦(») 


db)/2ir 


(2.3) 


is  the  temporal  coherence  function  of  the  light  field  and  the  Fourier  trans¬ 
form  of  the  spectral  density  4(b));  (0,  T)  is  the  Interval  during  which  elec¬ 
trons  are  counted.  We  assume  that  the  spectral  density  4(b))  Is  normalized 
so  that 


which  requires  that 


r-1 


f 


T^(0)  ■  T  p  4(u))  db)/  Ti  -  1. 


(2.4) 


(2.5) 
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In  terms  of  the  Fredholm  determinant 


D(u)  -  ][][  (1  + 
r“l 


(2.6) 


associated  with  the  Integral  equation  (2.2),  the  p.g.f.  Is 

00 

f(z)  -  [D(Np(l-z))r^  [l-X^Np(z-l)]-^  (2.7) 

r-1 

By  Introducing  the  residue  expansion  of  [D(u)l~^  we  can  write  this  p.g.f.  as 

f(2) 


r«l  P  r-1  '  ^  ' 

sfr 


N  \ 

\  “  T+nX**  ®*(“)  *  dD/du. 
P  r 


(2.8) 


Because  (1.23)  is  linear  In  the  derivatives  fj^,  the  probabilities  p^  can 
be  considered  as  a  weighted  s\im  of  probabilities  arising  from  the  individual 
terms  of  (2.8), 


p..  “  ✓  ®  p  » 

n  r  r,n* 


(2.9) 


r-1 


with 


n 


k«l 


5r- 


-G 


1  -  V 


l-y.e'®’  1-v  .-® 


(2.10) 
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Writing 


P 

r,n 


^n.k  “  S(n.  (2.11) 

and  using  the  recurrent  relation  (1.20)  for  the  modified  Stirling  numbers,  we 
find  a  recurrent  relation  for  the  terms  of  (2.11), 


(r) 


X  (e  T  +'r 

n+l,k  n+1  n,k-l  n,k  ^  * 


“  GCj.. 


1.1 


(2.12) 


The  sum  in  (2.9)  is  stopped  when  the  terms  decrease  to  the  point  of  insigni¬ 
ficance. 


(b)  Lorentz  Spectral  Density 

For  incident  light  with  a  Lorentz  spectral  density,  whose  normalized  form 
is 


«((») 


-  2m  , 
T(w^  +  p^) 


(2.13) 


the  p.g.f.  f(z)  of  the  primary  distribution  is  given  by  [17J 


f(z)  ■  e®Icosh w+§  ^ jsinh w] 

2  2 

w  -  m  +2Npm(l-z),  m  »  DT.  (2.14) 

We  denote  the  square  bracket  by  M(z)  and  write  it  in  terms  of  modified  spher¬ 
ical  Bessel  functions  ijj(w)  obeying  the  recurrent  relation 

‘n+lW  -  ^  l,(w)  <2.15) 

with 

i.l<w)  -  lo(»)  - 

whereupon 

1  2 
M(z)  -  2  ^-1^'')  +  Is 


r 


Then  following  Bedard  [18]  we  use 


/  N  m  r  . 

_  v 

2  -1 

(-  -S-j  [f  ♦ »(' 

-  (-1)\(2), 

which  defines  M  (2),  In  Leibnitz's  formula 


(2.17) 


_d 

dz 


r-O 


to  determine  a  recurrent  relation  for  the  derivatives  f  of  f(z). 

n  ’ 

fo  -  S  (r)  > 


r*0 


(2.18) 


fjj  “  f(z),  w  «  [m^  +  2Npm(l-2)l^. 


(2.19) 


The  probability  p^  of  counting  n  secondary  electrons  during  (0,  T)  when 
the  incident  light  has  a  Lorentz  spectral  density  is  then  calculated  by  sub¬ 
stituting  from  (2.19)  into  (1.23),  Alternatively,  one  can  calculate  the  prob- 

(1.17)  and  (1.2A)  by  Bedard's  method  [18],  which  corresponds 
to  taking  z  -  0  in  (2.14-.19), 


(c)  Negative  Binomial  Primary  Distribution 

The  primary  electron  distribution  can  often  be  closely  approximated  by 
the  negative  binomial  distribution 


n. 


(M) 

“ici 


^  (l-v)V, 


N 


N  +M’ 
P 


0,  1,  2,  . . . , 


(M)  -  M(M-»-l)...(M  +  k-l)  -  , 


'k  - - - 

In  which  the  number  M  of  degrees  of  freedom  is  given  by 


(2.20) 


11 


^•T\HO)rjj  (T-  |T|)I(^(T)rdT 


(2.21) 


in  terms  of  the  temporal  coherence  function  of  the  Incident  light  [16,  19] 

For  the  Lorentz  spectral  density  in  (2.13) 

<Kt)  -  (2.22) 


and 

M  ®  2m^/(2m- l  +  e**^®),  m  «  yT; 
M  1  as  m  -►  0,  M  -►  m  as  m  For  M  »  1, 


M  i  tU(0)1 


(2.23) 


(2.24) 


and  M  roughly  equals  the  time-bandwidth  product  WT,  where  W  is  the  equivalent 
bandwidth 


of  the  light.  The  p.g.f.  of  this  negative  binomial  distribution  is 


(2.25) 


f(2) 


N)’ 


and  the  derivatives  of  the  p.g.f.  f(z)  in  (1.23)  are 


'k  ■  <«k 


(2.26) 


(2.27) 


Again  using  the  recurrent  relation  (1.20)  with  (1.23)  we  can  write  the  prob¬ 
ability  p^  of  n  secondaries  as 


Pn"  PO 


k»  Pq  “  ^0’ 


(2.28) 


k«l 


2 


In  which  the  coefficients  T  .  obey  the  recurrent  relation 

Zl  ^  1C 

Vl.k  ■  TTT 


ve 


“G 


1  -  ve 


-G* 


»i,i  ■ 


(2.29) 


(d)  The  Neyman  Type-A  Distribution 

When  we  pass  to  the  limit  M  -*•  »,  keeping  the  mean  number  of  primary  elec¬ 
trons  fixed y  the  distribution  of  primary  electrons  turns  into  the  Poisson  dis¬ 
tribution, 

\  exp  (-Np)/k! ,  (2.30) 

for  which  the  p.g.f.  is 

f(2)  “  exp  tN  (z-1)];  (2.31) 

P 


and  the  distribution  of  secondaries  becomes  by  (1.18)  the  Neyman  Type-A  dls-  . 
trlbutlon  [2],  whose  p.g.f.  is 


h(z)  -  exp  [N  (e®^*“^^-l)]  (2.32) 

P 

and  for  which  the  probabilities  are  most  simply  computed  by  Neyman 's  recurrent 
relation  [20] 


^’k+l 


N  Ge  ® 

_E _ 

k  +  1 


Pk-r’ 


(2.33) 


r»0 

or,  for  n  »  1,  by  (1.24-.25)  and  (2.30). 

In  order  to  illustrate  the  transition  to  the  Neyman  Type-A  distribution 
with  Increasing  values  of  M,  ve  have  plotted  the  probabilities  p^,  n  >  0,  in 
Figs.  1  and  2  for  both  the  negative-binomial  primary  distribution  (solid  lines) 
and  the  distribution  arising  from  light  with  a  Lorentz  spectral  density  (dashed 
lines),  the  equivalent  value  of  m  having  been  determined  by  solving  (2.23). 


13 


Where  both  distributions  are  not  shown,  their  graphs  fell  too  close  to  be  dis¬ 
tinguished.  The  curves  marked  represent  the  Neyman  Type-A  distribution. 

In  Fig.  1,  *  2,  G  6;  in  Fig.  2,  >  6,  G  ■  2.  The  long  tails  of  the 

negative-binomial  distribution  for  small  values  of  M  carry  over  to  the  distri¬ 
bution  of  the  number  n  of  secondary  electrons.  The  negative-binomial  distribu¬ 
tion  is  seen,  furthermore,  to  yield  a  close  approximation  to  the  distribution 
.  arising  from  light  with  a  Lorentz  spectral  density  over  the  entire  range 

1  S  M  <  *,  the  t%j0  coinciding  at  the  extremes  M  =  1  and  M  =  except  that  the 
latter  drops  off  to  zero  slightly  less  rapidly  than  the  former. 


The  probabilities  p^  having  been  computed  by  these  methods,  the  cumulative 
probability  is  computed  by  the  summation  in  (1.5),  and  its  complement 
Is  determined  from  I-qJ*.  When  the  numbers  n  are  large,  however,  is  close 
to  1,  and  round-off  errors  corrupting  the  probabilities  p^  Introduce  large  rela- 
tlve  errors  into  the  tall  probability  .  These  methods  furthermore  require 
storage  of  more  and  more  numbers  and  entail  more  and  more  additions  and  multipli¬ 


cations  as  n  and  E(n)  Increase.  We  therefore  turn  to  methods  that  enable  computa¬ 
tion  of  the  tall  probabilities  and  directly  and  are  the  more  efficacious, 


the  larger  the  mean  number  N 


of  primaries  and  the  gain  G  in  the  multiplicative 


process . 
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III.  The  Method  of  Residues 


for  all  integers  k,  -»  <  k  <  «*, 

As  shown  in  Appendix  A,  we  can  expand  the  contour  C  into  a  rectangle  at 
infinity,  provided  that  we  also  enclose  each  pole  hy  a  small  circle  tra-* 

versed  clockwise.  The  integral  around  the  rectangle  vanishes,  and  we  are  left 
with  the  results  of  applying  the  residue  theorem  to  the  integrals  around  each 
pole,  whereupon 


(3.5) 
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It  Is  furthermore  shown  In  Appendix  A  that  the  error  Incurred  by  stopping  the 
suanation  in  (3.5)  at  k  «  K  is  bounded  by 

1Rj,|  «  ir  (1-v^)Yjj^^^2k  +  1)]  • 

\.  :  v  -  Tr^r(S^/r(n/2).  (3.7) 


The  factors  remain  bounded  and  are  of  order  n”^  for  n  »  1.  When  n  is  large, 
only  a  few  terms  contribute  significantly  to  the  sum  in  (3.5),  and  convergence 
Is  rapid. 

For  the  Lorentz  spectral  density  in  (2.13),  to  which  we  have  applied  this 
method,  the  eigenvalues  are  given  by 

\  2*  ® 

+m 

where  the  c^  are  the  solutions  of  the  transcendental  equation 

2  2 

cot  ”®  )/2mc^.  r  ■  1,  2,  ..., 

which  can  also  be  written 

c^  tan  (c^/2)  ■  m,  r  odd 

Cj,  cot  (c^/2)  •«  -m,  r  even  (3.9) 

The  quantity  c^  lies  between  (r“l)ir  and  nr,  and  with  c^  *  (r-l)7r  +  E^,  we  find 
for  r  »  m/ir. 


s:  2m/rTi. 


(3.10) 


The  coefficients  in  (2.8),  (2.9),  and  (3.1)  are  given  by  [21] 

r  .  .  r\  sin  c^ •  r  I  ’ 


which  after  some  algebra  can  be  reduced  to 

o,  -  (-l)'"^e®X^(2-mX  )/(l  +  X  ). 


(3.11) 


For  large  values  of  r,  by  (3.11)  and  (3.8) 

|Oj.|  s-  4oe“/r^ir^,  (3.12) 

_2 

which  decreases  in  proportion  to  r  as  r  This  rapid  decrease  does  not 

set  In,  however,  until  r  exceeds  m/TT.  When  m  »  1,  the  coefficients  for 
r  <  n/ir  are  large  In  absolute  value;  they  always  alternate  In  sign.  When  the 
number  n  Is  of  the  order  of  E(n)  or  larger,  the  convergence  of  (3.1)  Is  acceler- 
ated  by  the  rapid  decrease  of  the  terms  In  (3.5)  with  Increasing  r,  but 

for  n  «  E(n)  many  terms  of  the  series  must  be  taken  when  m  »  1.  Because  the 
terms  of  (3.1)  are  of  decreasing  magnitude  and  alternate  In  sign;  the  error  is 
bounded  by  the  last  term  Included  In  the  sum. 

In  Table  I  we  compare  the  results  of  using  (3.1)  and  (3.5)  with  those  ob¬ 
tained  by  summing  the  probabilities  calculated  by  (1.23)  with  (2.19)  in  double 
precision.  The  column  headed  "last  Increment"  lists  the  last  term  added  Into 
the  sum  In  (3.1).  We  took  m  «  4.4365821,  corresponding  to  M  -  5  in  (2.23),  and 
forty  terms  of  (3.1)  were  summed.  For  these  we  found 


40 

Q*  -y^  -  1.046101 

r-1 

Instead  of  1.  It  Is  seen  that  the  error  Incurred  by  the  residue  series  decreases 
with  Increasing  n,  and  the  faster,  the  smaller  the  gain  G. 

When  the  p.g.f.  f(z)  of  the  primary  electrons  can  be  approximated  as  in 
(2.26)  for  integral  values  of  M,  the  p.g.f. 

possesses  a  vertical  row  of  poles  of  order  M  at  the  points 

C|^  -  l  +  G“^(lnv“^  +  2kTrl),  -«  <  k  <  •.  (3.1A) 
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Expanding  the  contour  of  Integration  across  this  row  of  poles  and  applying  the 
residue  theorem  as  before  yields  for  the  complementary  cumulative  probability 


n+  -1  d**~^  r.  _  .M  s~”h(z)1| 

It  Is  shown  in  Appendix  B  that  this  can  be  written  as 

-n  M-1 


(3.15) 


(1-v) 


’n'  -  E 


k— a 


j*0 


kj 


(3.16) 


with 


(")- 


'oil. 

/  ^  si 
S“0 


-j 'on. 


n(n  +  l)...(n  +  s-  1), 


(3.17) 


(3.18) 


(M) 


and  the  coefficients  are  tabulated  In  Table  II.  The  manner  of  calculating 

these  coefficients  Is  given  In  Appendix  B.  Again  the  symmetry  of  the  poles  about 
the  real  axis  permits  carrying  the  summation  from  k  «  0  to  and  replacing 
terms  with  k  >  0  by  twice  their  real  parts. 

Table  III  compares  the  results  of  this  residue  series  with  the  exact  prob¬ 
abilities  for  H  ■  5  and  the  same  parameters  as  in  Table  I.  The  column  labeled 
'Hum*  lists  the  maximum  value  of  k  In  eq.  (3.16),  In  which  the  summation  was 
stopped  when  the  ratio  of  the  absolute  value  of  the  last  term  to  that  of  the 
sum  fell  below  10  This  Is  of  the  order  of  smgnltude  of  the  discrepancy  be¬ 
tween  the  exact  values,  computed  by  (2. 28 -.29),  and  those  computed  by  (3.16). 


As  is  to  be  expected  from  (3.6),  the  number  of  terms  to  be  summed  decreases  with 
Increasing  n.  The  column  of  Table  III  headed  "Single  term"  lists  the  contribution 
of  the  pole  Cq  on  the  real  axis  alone:  for  n  greater  than  E(n)  It  Is  seen  to 
provide  an  accurate  approximation  to  the  right-hand  tail  probability  Q^. 


IV.  Numerical  Contour  Integration 


Ihe  method  of  residues  In  Sec.  Ill  cannot  be  applied  to  the  Neyman  lype-A 

distribution,  whose  p.g.f.  has  no  poles,  nor  to  the  distribution  of  secondaries 

arising  from  primaries  with  a  negative-binomial  distribution  for  a  nonintegral 

value  of  M,  for  which  the  singularities  are  branch  points.  For  light  with  a 

spectral  density  such  as  the  Lorentz  in  (2.13),  furthermore,  (3.1)  requires  a 

great  many  terms  when  WT  »  1  and  n  «  E(n)  In  order  to  determine  Q  ”  *  1  -  Q 

n  n 

with  usefully  high  relative  accuracy.  We  therefore  resort  to  computing  the  tail 
probabilities  and  by  evaluating  the  contour  Integrals  in  (1.9)  and  (1.10) 
numerically. 


For  the  sake  of  efficiency  one  would  like  to  Integrate  along  that  contour 
C  on  idilch  the  magnitude  of  the  Integrand  decreases  as  rapidly  as  possible  from 
Its  maximum  value,  which  occurs  for  z  real  and  positive.  Such  a  path  Is  known 
as  the  path  of  steepest  descent  [22].  With  the  Integrand  written  In  the  form 

exp['l'(z)l  «  ±z‘’"h(z)/(z  -  1),  (4.1) 

the  imaginary  part  Im'P(z)  of  the  "phase”  V(z)  Is  constant  along  this  path.  The 
path  of  steepest  descent  furthermore  passes  through  the  saddlepolnts  of  the  In¬ 
tegrand,  which  are  the  points  at  which 


dz 


In  h(z) 


n  _  1 

z  ”  z  -  1 


0, 


(4.2) 


If  one  plots  the  magnitude  of  the  Integrand  exp['i'(x)l  for  values  of  z  x 


on  the  positive  real  axis,  one  finds  that  it  Is-  a  convex  function  with  one  mln- 
isnim  at  a  point  Xq  in  0  <  x  <  1  and  another  at  a  point  x^  in  1  <  x  <.  Y^,  where 
Yq  la  the  leftmost  singularity  of  h(z)  on  the  real  axis.  These  points  z  «  x^ 
are  roots  of  (4.2)  and  the  principal  saddlepolnts  of  the  Integrand.  When  the 
contour  of  integration  passes  vertically  through  z  xj*  or  z  ■  x^  ,  the  magni¬ 
tude  of  the  integrand,  maximum  at  the  saddlepolnt,  decreases  most  rapidly  on 
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either  side. 


Figure  3  exhibits  typical  paths  of  steepest  descent  for. (1.9)  and  (1.10) 
vhen  these  are  used  to  calculate  the  Neyman  Type-A  distribution,  whose  p.g.f. 
h(*)  Is  given  In  (2.32).  They  are  drawn  for  •  G  ■  10.  Only  the  curves  In 
the  upper  half-plane  are  Illustrated;  the  portions  In  the  lower  half-plane  are 
their  Blrror  linages.  The  left-hand  set  of  curves  refers  to  for  n  =  75,  the 
right-hand  set  to  for  n  »  150.  Small  circles  Indicate  the  saddlepolnts. 

The  paths  of  steepest  descent  go  off  to  Infinity  along  asymptotes  at  values  of 
y  -  lo  z  equal  to  odd  multiples  of  tt/g. 

Utilizing  the  path  of  steepest  descent  would  require  computing  a  number  of 
saddlepolnts  Zq  by  solving  (4.2)  and  then  tracing  the  branch  of  the  path  of 
steepest  descent  passing  through  each,  a  cumbersome  procedure.  We  therefore 
instead  chose  as  our  contour  of  integration  a  vertical  straight  line  passing 
through  the  saddlepolnt  xj  for  and  through  x^  for  Q^.  This  line  passes 
close  to  the  saddlepolnts  of  the  Integrand  lying  above  and  below  the  real  axis. 
(The  advantages  of  Integrating  along  a  path  passing  through  or  near  a  string  of 
saddlepolnts  were  pointed  out  by  Lugannanl  and  Rice  [23].) 

The  saddlepolnts  x^  and  on  the  real  axis  can  most  expeditiously  be 
found  by  solving  (4.2)  by  Newton's  method,  starting  with  an  initial  trial  value 
just  to  the  left  of  z  •  1  for  x^  and  just  to  the  right  of  z  »  i  for  Xq  ;  at 
each  stage  one  replaces  the  trial  value  Xq  by 

x^  Xq’  -Y’(xq’ )/r'(xQ’ ),  (4.3) 

wlthY'  the  first,  T'  the  second  derivative  of  the  phase  V(x);  see  (4.1).  Along 
a  vertical  contour  through  xj  or  x^  the  magnitude  of  the  integrand  In  (1.9)  or 
(1.10)  decreases  most  rapidly.  For  n  <  E(n)  it  Is  most  expeditious  to  evaluate 
In  (1.9)  by  deforming  the  contour  Into  a  straight  line  through  the  left-hand 
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_  if 

saddlepolnt  Xq  ;  for  n  >  E(n)  one  evaluates  by  (1.10)  and  deforms  the  con¬ 
tour  Into  a  straight  line  through  the  right-hand  saddlepolnt  x^. 

For  reasons  discussed  In  [24]  the  trapezoidal  rule  Is  recommended  for  In¬ 
finite  Integrals  of  analytic  functions.  For  z  >  x^  +  ly  In  the  neighborhood  of 

the  saddlepolnt  Xq  »  x^  or  x^'  the  Integrand  has  approximately  a  Gaussian  de- 

-4 

pendence  on  y  with  a  width  of  the  order  of  [¥''(Xq)]  Using  the  trapezoidal 
rule  In  the  form 


where  -  (Xq^  +  Yo^)^* 

For  the  Neynan  Type-A  distribution, 

I  I  G(XfL-l) 

max  |h(x_  +  ly)l  •  exp  [N  (e  °  -1)] 

y  ^ 

end  for  that  arising  from  a  negative-binomial  primary  distribution,  by  (2.26) 
and  (1.18), 

max  |h(xQ  +  iy)|  -  (1  -  v)“ll  -  v  exp  G(Xq  -  1))"“. 

y 

The  final  bracket  Is  positive  because  the  saddlepoint  x^  lies  to  the  left  of  the 
leftmost  pole  of  h(z)  as  given  by  (3.14)  %rlth  k  -■  0. 

Tables  IV  and  V  show  results  of  our  computation  of  the  Neyman  Type-A  distri¬ 
bution  and  of  that  arising  from  a  negative-binomial  primary  distribution  with 
M  "  5.  The  summations  were  halted  when  the  ratio  of  the  bound  on  the  truncation 
error,  as  given  by  (4.5),  to  the  computed  probability  fell  below  10  The 
column  headed  "Exact”  was  computed  by  summing  probabilities  p^  calculated  In 
double  precision  from  (1.17)  with  (2.30)  and  (2.20)  respectively,  and  because 
of  round-off  error  In  our  computer,  which  even  In  double  precision  carries  only 
about  sixteen  decimal  digits,  the  nixmerlcal  Integrations  had  to  be  carried  out 
in  double  precision  as  well.  The  numerical  contour  Integration  in  these  examples 
yielded  the  tail  probabilities  to  six  significant  figures  with  fewer  than  two 
hundred  steps. 

In  Figs.  4  and  5  we  exhibit  the  cumulative  distribution  for  the  parameters 

used  in  Tables  IV  and  V:  N  ■5,G«20  and  N  •°20,  G-5,  respectively,  and 

P  P 

for  M  ■*  1,  3,  5,  10,  and  *,  the  mean  N^G  *  100  remaining  fixed.  The  figures  il¬ 
lustrate  the  manner  In  which  the 'Neyman  Type-A  distribution  (H  >  *»)  is  approached 
with  increasing  M.  Figs.  6  and  7  exhibit  the  same  cumulative  distributions,  but 

for  H  ■  18,  G  ■  72  and  -  72,  G  ■  18,  respectively,  rhe  mean  N  G  -  1296 

P  P  P 
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Tenalnlng' fixed.  Comparison  with  Figs.  A  and  5  shows  the  approach  of  the  Neyman 
lype-A  distribution  to  normality. as  the  mean  N^G  Increases,  as  predicted  by  Telch 
[2],  The  distributions  arising  from  the  negative-binomial  primary  distribution 
manifest  no  such  progression  toward  normality. 

The  small  circles  in  Figs.  4  to  7  mark  the  values  of  the  cumulative  distri¬ 
bution  of  secondary  electrons  arising  from  incident  light  with  a  Lorentz  spectral 
density:  these  were  computed  by  the  residue  series  in  Sec.  III.  For  a  few  sets 
of  values  the  numerical  integration  method  of  this  section  was  also  carried  out, 
and  results  agreed.  The  tails  of  this  distribution  do  not  drop  off  to  zero  so 
rapidly  as  for  the  approximating  distribution  calculated  from  the  negative-binomial 
primary  distribution  with  the  same  number  M  of  degrees  of  freedom,  as  specified 
by  (2.23).  The  latter  corresponds  to  a  spectral  density  that  cuts  off  sharply 
at  a  frequency  deviation  of  the  order  of  W  =  M/T  from  the  central  frequency  of 
the  light;  the  Lorentz  spectrum,  on  the  other  hand,  has  very  long  tails.  Because 
very  large  frequency  deviations  are  much  more  prevalent  in  the  latter,  there  are 
more  opportunities  for  large  count  deviations  to  be  accumulated,  and  numbers  n 
much  larger  or  much  smaller  than  the  mean  E(n)  are  therefore  more  probable  than 
with  primaries  having  a  negative-binomial  distribution. 


V.  Saddlepolnt  Approximations 


The  principal  contributions  to  the  Integrals  in  (1.9)  and  (1.10),  after 
the  contours  have  been  deformed  into  the  paths  of  steepest  descent  as  discussed 
in  Sec.  IV,  come  from  the  neighborhood  of  each  saddlepolnt  these  saddlepoints 
are  the  roots  of  (4.2).  The  Integrals  can  then  be  approximated  by 


00 

0^=^  a  +2  Re^  [2Trr'(Zj^)r^exp  [S' (z^)  ]  (1  +  C^)  ,  (5.1) 


k=l 


±(0) 


tidiere  represents  the  contribution  of  the  branch  of  the  path  of  steepest 

descent  passing  through  a  saddlepolnt  on  the  real  axis;  z^,  Z2>  saddle- 

points  lying  above  the  real  axis,  and  is  a  complex  correction  term  given  by 


’l'^“^(z,j)/tY^^^(Zv)l“'^^,  «i'^"*^(z)  -  d“4'(2)/d2“. 


(5.2) 


In  (5.1)  we  must  take  Re  [^"(z.  )]  >  0.  Since  the  contributions  of  the  complex 

±(0) 


saddlepoints  are  usually  much  smaller  than  it  is  often  unnecessary  to 

Include  the  correction  Cj^.  Occasionally,  however,  Cj^  may  have  a  strong  influence 


if  the  principal  term  that  (1  +  Cj^)  multiplies  has  a  phase  close  to  ±it/2,  so  that 


without  the  factor  (1  +  C^)  its  real  part  would  be  small.  The  saddlepoints 
must  be  determined  with  high  precision. 


One  method  of  approximating  the  on-axls  term  Q^' 
part  of  (5.1-.2)  [9], 


±(0) 


is  to  use  the  counter- 


[27rf^2>(xQ)]“^expf(xQ) 


(5.3) 


with  X. 


x,r  for  Q  and  x. 


+(0) 


for  Q  ,  where  x.  and  x^  are  the 


real*-axis  saddlepolnts  defined  In  Sec.  IV.  This  approximation  Is  the  more  ac¬ 
curate,  the  farther  the  value  of  n  lies  from  the  mean  E(n)  of  the  distribution. 

For  n  near  E(n)  greater  accuracy  is  achieved  by  utilizing  the  first  two 
terms  of  a  uniform  asymptotic  expansion  [12,  25].  It  is  based  on  the  modified 


phase 


y(z)  =  In  h(z)-nln2. 


of  which  a  saddlepolnt  x^  is  the  root  of 


t’Cz) 


h'(2)  n  ^  n 
h(2)  "  z  ^ 


(5.4) 


This  saddlepolnt  z  =  Xq  can  also  be  quickly  computed  by  Newton's  method.  There 
is  a  single  such  saddlepolnt  on  the  real  axis,  0  <  x^  <  Yq»  where  Yq  Is  the  left- 
most  singularity  of  h(z)  on  the  real  axis*  At  the  point  x  »  'i'(x)  Is  minimum. 

The  contour  of  integration  is  displaced  from  the  path  of  steepest  descent 
of  the  entire  Integrand  so  that  It  coincides  with  that  of  the  function  exp'FCz), 
lAlch  lies  nearby.  Then  the  first  two  terms  of  the  uniform  asymptotic  expansion 
[25]  approximate  the  tail  probabilities  as 


n  >  E(n); 


'iui  /«-  \  1  B 


[-2'f(x0)]" 

n<E(n):  M 

+  {|xq-  ir^[2TT?"(7Q)]"'^-  i  [-7r?(xQ)]"^}  exp  [?(xq)],  (5.5) 


in  which  primes  indicate  differentiation  and 


erfc  u  ■  (2it)' 


exp  (-t^/2)  dt 


la  the  error-function  integral.  Corrections  of  higher  order  are  listed  in  [25]. 
For  n  E(n)  >  N^G,  Xq  ■■  1  and  this  approximation  breaks  down,  but  by  Eq.  (17) 
of  [25]  we  can  write  for  n  *  E(n)  ■  n 


(5.6) 


further  terms  are  to  be  found  in  [25],  Thus  for  the  Neyman  Type-A  distribution 

(5.7) 

and  for  the  distribution  of  secondaries  generated  by  primaries  with  a  negative- 
binomial  distribution 

.  i  -  i  j  .  V)  .  7,  j-^  . 

V  =  N  /M.  (5.8) 

P 

In  Tables  VI  and  VII  we  compare  the  results  of  the  saddlepolnt  approxima¬ 
tion  with  the  exact  values  of  the  cumulative  distribution;  Table  VI  refers  to 
the  Neyman  Type-A  distribution.  Table  VII  to  that  arising  from  the  negative- 
binomial  primary  distribution.  The  column  marked  "UAE"  was  computed  from 

(5.5)  or  —  for  n  =  N  G  --  from  (5,6).  The  "off-axis  contribution"  is  the  sum- 

P 

matlon  in  (5.1)  over  the  four  nearest  complex  saddlepoints  z^,  Z2»  Zy  and  z^, 
including  the  correction  which  amounted  to  a  few  percent.  This  contribution 
is  the  larger,  the  smaller  n;  for  n  >  E(n)  it  is  hardly  significant.  The  smaller 
the  gain  G,  the  farther  the  complex  saddlepoints  lie  from  the  real  axis,  and  the 
less  they  contribute  to  the  total  probability  or  Q^.  The  column  labeled 
"crude  SP"  lists  the  value  of  computed  from  (5.3),  omitting  the  factor  in 

brackets.  This  is  the  most  simply  calculated  of  approximations  to  these  prob¬ 
abilities  and  is  often  adequate  when  n  lies  far  In  the  tails  of  the  distribution 
N  »  1,  G  »  1.  It  can  be  useful,  for  Instance,  in  initial  search  for  the  level 
n  yielding  a  particular  value  of  or  Q^,  as  in  setting  a  bias  level  to  attain 

a  pre-asslgned  false-alarm  probability,  or  "size",  in  a  hypothesis  test. 
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Conclusion 


Three  methods  have  been  presented  for  calculating  the  cumulative  distribu¬ 
tion  0^  and  the  complementary  distribution  of  the  number  n  of  secondaries 
In  a  particle-multiplication  process.  They  are  appropriate  when  the  number  n  and 
its  mean  value  E(n)  »  N^G  are  so  large  that  the  recurrent  computations  of  Sec.  I 
are  Infeasible.  When  the  probability  generating  function  h(z)  has  poles  whose 
locations  and  residues  can  be  determined,  the  residue  series  of  Sec.  Ill  is  the 
most  efficient  for  n  ^E(n).  For  n<E(n)  in  these  cases,  and  in  general  when 
a  residue  expansion  cannot  be  applied  —  as  for  the  negative-binomial  primary  dis¬ 
tribution  with  M  nonintegral,  and  for  the  Weyman  Type-A  distribution  —  the  nu¬ 
merical  integration  method  of  Sec.  IV  can  be  made  as  accurate  as  desired  by  using 
sufficiently  small  steps.  If  approximations  suffice,  the  saddlepoint  methods  of 
Sec.  V  can  be  utilized,  but  one  must  in  general  include  the  contributions  of 
aaddlepolnts  above  and  below  the  real  axis  in  addition  to  that  of  the  principal 
saddlepoint  on  the  real  axis. 
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Appendix  A.  Convergence  of  Residue  Series 


We  expand  the  contour  In  (3.2)  Into  a  rectangle  with  vertical  sides  at 
X  >  and  with  horizontal  sides  along  the  lines  z  x ± iy  ,  where 

Iv 

-  2ir(K  +  l)/G, 

<  X  <  ®.  (A.l) 

Cr^ 

In  the  course  of  its  expansion  the  contour  crosses  the  line  of  poles  at  , 

-K  ^  k  ^  K,and  leaves  behind  It  a  little  circle  surrounding  each  pole. 

The  vertical  sides  of  the  rectangle  contribute  zero  to  the  Integral  around 
it  because  of  the  factor  z  "  In  the  Integrand.  The  contributions  from  the  little 
circles  are  evaluated  by  the  residue  theorem  as  in  Sec.  3.  and  we  obtain  for 


(3.2) 


trlth 


where 


q  -  q  +IV 

^r,n  ^r,n  TC 


1-v  K  r(r) 

(K)  .  r  ^  ^ 

*r  ,n  G 


-n 


I _ V  * 


k=-K  ^k 


Z-1  27Ti« 


-r 

•'®fiy^ 


iz_ 

l  -  1  2111  ■  V  • 

We  want  to  bound  the  remainder  term  R|^. 

Along  the  upper  side  of  the  rectangle,  z  ■  x  +  ly^^,  -»  <  x  <  ®, 

^G(z-l)  _  exp[G(x-l) +2TT(K  +  i)i]  - 


(A.  2) 


(A.  3) 

(A.4) 


(A.  5) 


(A.6) 


and  since 


we  find 


1^  2Tr(x- l  +  lyj^)(x  +  lyj^) 

^  q-v^)  r _ ^ _ 

la,  U-i  +  iyKl(*^+yK^)"^ 


|x-l  +  lyj^r^  S 


1-v  f* 


.  2^  2.-n/2  . 

(x  +yj^  )  dx 


h''-  -"  f 

2i.  / 


cos”"^  0  d0 


(A.  8) 


-  “IF-  Vk 

with  given  by  (3.7).  Thus  the  remainder  term  Is  bounded  by 


-  |2RelJ  S  2|lJ 


w  ^<1-Vj.)y^^^2k  +  1))  • 


(A.  9) 


(A. 10) 


The  same  bound  applies  to  the  error  Incurred  by  cutting  off  the  summation 

M 

in  (3.16)  at  k  >  -K  and  k  ■  K,  with  (1-v^)  replaced  by  (1-v)  ,  as  can  easily 
be  shown  by  replacing  h^(z)  In  (A. 7)  with  h(z)  In  (3.13). 

In  this  way  we  show  that  the  partial  sums  In  (A. 3)  converge  to  the  prob¬ 
abilities  q  _  defined  In  (3.2)  and  summed  In  (3.1).  The  convergence  of  the 


series  In  (3.1)  follows  from  Abel's  criterion  [26].  The  series 

00 

r-1 

converges  to  1  by  virtue  of  (2.8)  with  z  -  1;  f(l)  -  1  by  (2.7).  We  need  only 
to  show  that  for  r  sufficiently  large,  the  ^  form  a  monotonely  decreasing 
sequence  In  r  for  fixed  n.  Now  by  (2.20),  (2.26),  and  (1.17)  with  v  -  v^, 

M  «  1,  the  quantities  P  defined  by  (2.8)  and  (2.9)  are 

A 

00 

*‘r,m  *  -  ^r>  ^  v^*^(kG)”  (A.12) 

k-1 

The  eigenvalues  of  (2.2)  and  hence  the  v^  defined  In  (2.8)  are  arranged  in 
descending  order,  and  v^  0  as  r  ^  For  some  integer  r^,  Vj.<  i  for  all 
r  >  r^,  whereupon 

(1  -  >  (1  -  I  >  r„,  *k  »  I, 

and  the  probabilities  (l**v^)v^  In  (A.12)  form  a  monotonely  decreasing  sequence. 
The  probabilities  therefore,  form  such  a  sequence,  and  so  do  the  prob¬ 

abilities  ' 

«» 

q  -Tp 

^r,n  r,m 

k*m 

Since  the  series  in  (A. 11)  converges,  and  for  r  >  r_  the  q  form  a  monotonely 
decreasing  sequence,  the  series  In  (3.1)  converges  by  Abel's  criterion. 
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Appendix  B.  Residue  Expansion 


Dropping  the  subscript  k,  we  write  the  terms  In  (3. 15),. with  (3.13),  as 
Res 


with 


so  that 


with 


(M 


-1  r  11 

^  ^  [r(.)  c(o]| 


Hr 

(M 


F(*) 


G(z) 


-M 

_2 _ 

Z  “  1 


Res 


Hr 

(M 


_ _ f\  ttZ  oz 


r*0 

M-1 


^1-r 

dz^ 


z=C 


r«0 


♦r  -  ^  77 


^r  -  ^  77 


z-5 


Z-t 


Starting  with  we  put  Into  F(z)  In  (B.2) 

G(z-C)  -  y, 

80  that 


(B.l) 


(B.2) 


(B.3) 


(B.4) 


(B.5) 


(B.6) 
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32 


\  ’rrr  \ — 7*  l'^^:’■  .*-'  .*»  ■:■  •?  r^r^^r^r^r-'.' 


H  -_i-  V'  ”m.k 

r+l.k+l  r  +  1  ^  ir  -  m) !  (r  +  2  -  m) 


B«k 

r-k 


*  **  U 

^  ”r 

+  1  jl(, 


r  +  1  jI(j+2) 

j-0 


with  special  values 


"l.l  ■  r  “t.l  ■  (r+\)!-  ®r.r  '  TT: 

2  rf 


For  In  (B.5)  we  obtain 


Tr-^  77  1)“" 

dz 


z-C 


s*0 


T-8 
dz* 


-  c-i)’'  V 


It  Is  convenient  to  write 

-»  .  (-l)“-l-V^G-<"-'-«I<« 
r  r  * 

I «)  .  (-i)r  ^ 


with 


k-1 


(M) 


1,  0  ^  r  s  M-  1, 


whereupon,  with 
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(B.13) 


(B.14) 


(B.15) 


(B.16) 


(B.17) 


III  tim  »  il  ■  .•|i<..i<l  .'•i  'm,mtim  ifc.  li.  J!i' .  ■• '. 


(B.18) 


ve  find  from  (B.3) 


8«0 


Res 


r«0 


(1  -  V)**  ^  (M)  -j 

G~  2^°j  C-1 

j-0 


E. 


a-D-J  V 

8«0 


1.  («). 


-j  '“'8  X,  --ivS 


s! 


(1-C 


where  the  coefficients 


-  (M) 


P  S  j  i  M-1 


(B.19) 


(B.20) 


(B.21) 


are  tabulated  In  Table  II.  A  residue  of  this  form  Is  calculated  for  each  pole 
5  ■  and  the  result  Is  the  ctxpresslon  (3.16)  for  the  probability  Q^. 
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Table  I 


Performance  of  Residue  Series  for  (40  terms) 
Lorentz  spectral  density,  m  4.4365821 


■ 

N  - 
P 

5.  G  » 

20 

N  > 
P 

20,  G  » 

5 

n 

Residue 

series 

Exact 

Last 

Increment 

n 

Residue 

series 

Exact 

Last 

increment 

20 

0.940032 

0.939034 

2.20(-3) 

20 

0.991725 

0.991717 

6. 38 (-6) 

0.842027 

0.842012 

2. 40 (-5) 

40 

0.928487 

0.928483 

3. 60 (-10) 

0.710186 

0.710182 

3.45(-7) 

60 

0.785282 

0.785280 

2. 03 (-14) 

80 

0.567215 

0.567213 

6. 71 (-9) 

80 

0.602037 

0.602037 

1.14 (-18) 

00 

0.433256 

0.433255 

1.25 (-10) 

100 

0.428268 

0.428268 

6.43(-23) 
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Table  II 


Coefficients 

in  (3.16) 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

1 

2 

1 

35 

5 

1 

12 

2 

X 

15 

4 

17 

4 

3 

1 

203 

49 

35 

7 

1 

45 

8 

16 

2 

1 

469 

967 

28 

23 

A 

90 

120 

3 

3 

X 

29531 

801 

1069 

27 

39 

9 

1 

5040 

80 

80 

2 

4 

2 

X 

6515 

4523 

285 

3013 

75 

145 

e 

1008 

378 

16 

144 

4 

12 

X 

177133 

84095 

341693 

8591 

7513 

605 

44 

11 

1 

25200 

6048 

1512 

288 

240 

24 

3 

2 

X 

190553 

103344293 

1254429 

242537 

1901 

10831 

35 

A 

1 

25200 

6531840 

45360 

6048 

40 

240 

«X3 

2 

Table  III 


Performance  of  Residue  Series  (3.16)  for  Q^'*’ 
Negative-Binomial  Primary  Distribution  (M  =  5) 

Single  term  Residue  series  Num  Exact 

Np  •=  5.  G  -  20 


20 

0.932417 

0.955457 

4 

0.955457 

40 

0.833080 

0.873738 

3 

0.873737 

60 

0.705184 

0.749536 

2 

0.749535 

80 

0.568492 

0.601921 

2 

0.601920 

100 

0.439202 

0.439582 

2 

0.439582 

150 

0.198436 

0.198413 

2 

0.198413 

200 

0.0763843 

0.0763844 

2 

0.0763845 

N  = 
P 

20,  G  »  5 

20 

0.983799 

0.983798 

1 

0.983799 

40 

0.911155 

0.911153 

1 

0.911155 

60 

0.775914 

0.775912 

1 

0.775914 

80 

0.608256 

0.608256 

1 

0.608258 

100 

0.443583 

0.443581 

1 

0.443583 

150 

0.157956 

0.157956 

1 

0.157956 

200 

0.0440489 

0.0440489 

1 

0.0440489 
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Table  IV 


;  Neyman  Type-A  Distribution 
Numerical  Contour  Integration 


Number  of 

^0 

Result  of 

steps 

u 

Integration 

N  - 

5,  G  «  20 

P 

37 

1.687 

3.5527190(-2) 

74 

1.687 

3.5523618(-2) 

148 

1.687 

3.5523618(-2) 

36 

0.6551 

0.52720529 

71 

0.6460 

0.52698250 

141 

0.6415 

0.52698245 

36 

0.5099 

1-6. 3492573 (-2) 

71 

0.5029 

l-6.3492302(-2) 

141 

0.4993 

l-6.3492302(-2) 

N  - 
P 

20,  G  -  5 

41 

2.856 

9.3933963 (-5) 

82 

2.856 

9.3933799(-5) 

164 

2.856 

9.3933799(-5) 

30 

0.9178 

0.51061370 

60 

0.9178 

0.51045010 

120 

0.9178 

0.51045007 

25 

0.7618 

1-2. 9833051 (-3) 

50 

0.7618 

1-2. 9833045 (-3) 

99 

0.7542 

l-2.9833045(-3) 

Exact 

3.5523618(-2) 

0.52698245 

1-6. 34 92302 (-2) 

9.3933799(-5) 

0.51045007 

1-2. 9833046 (-3) 


Table  V 


Distribution  Generated  by  Negatlve^Blnomlal,  M  »  5 


Numerical  Contour  Integration 


n  Number  of  y. 

steps 


Result  of 
integration 


Exact 


Np  •=  5,  G  =  20 


25 

30 

1.3191 

1-9. 7678911 (-2) 

60 

1.3191 

1-9. 7652467 (-2) 

119 

1.3081 

1-9. 7652467 (-2) 

100 

71 

0.5602 

.0.43963053 

143 

0.5563 

0.43958206 

285 

0.5543  • 

0.43958206 

175 

57 

0.4193 

0.12519049 

115 

0.4157 

0.12518192 

229 

0.4157 

0.12518192 

Np  -  20,  G  »  5 


25 

22 

1.3094 

1-2. 7957485 (-2) 

43 

1.2796 

1-2.7954584 (-2) 

84 

1.2647 

l-2.7954585(-2) 

100 

55 

0.5571 

0.44363275 

112 

0.5571 

0.44358254 

223 

0.5571 

0.44358254 

175 

47 

0.4211 

8. 5479092 (-2) 

94 

0.4211 

8.5472726(-2) 

187 

0.4188 

8.5472726(-2) 

l-9.7652465(-2) 


0.43958206 


0.12518192 


1-2. 7954583 (-2) 


0.44358255 


8.5472739(-2) 


Table  VI 


SaddlepolnC  Approximations 


Neyman  Type-A 


N  -  5,  G  -  20 
P 


n 

UAE 

Off-axis 

contribution 

Total 

Exact 

25 

3.14560 (-2) 

4.09137 (-3) 

3.55474(-2) 

3. 55236 (-2) 

50 

0.130576 

9. 51063 (-6) 

0.130585 

0.130600 

100 

0.527499 

-5. 51718 (-4) 

0.526947 

0.526982 

150 

0.857731 

3. 01480 (-5) 

0.857762 

0.857777 

200 

1-2. 50472 (-2) 

'  9. 23185 (-7) 

1-2. 50463 (-2) 

1-2. 50431 (-2) 

N  -=  18 
P 

,  G  »  72 

250 

1.98051(-5) . 

-1.71780 (-6) 

1.80873(-5) 

1.80638 (-5) 

500 

1.8l082(-3) 

1.81503 (-5) 

1.82897 (-3) 

1.82849(-3) 

750 

2. 98581 (-2) 

5. 45979 (-5) 

2. 99127 (-2) 

2. 99100 (-2) 

1000 

0.168343 

-2. 72697 (-4) 

0.168070 

0.168064 

1500 

0.752654 

-8.41239(-5) 

0.752570 

0.752567 

2000 

1-1. 58622 (-2) 

-2. 85533 (-7) 

1-1. 58625 (-2) 

1-1. 58627 (-2) 
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•**  '*'^'**'*  •  •  ^  »,•  •, 


Crude  S-P 

3. 289 (-2) 
0.1324 
0.5103 
0.8626 
l-2.478(-2) 


2. 022 (-5) 
1.826 (-3) 
2. 986 (-2) 
0.1657 
0.7648 
1-1. 571 (-2) 


Table  VII 


Saddlepolnt  approximations 

:  Distribution  Arising  from  Negative-Binomial,  M  =  5 

N  -  5.  G  -  20 
P 


n 

UAE 

Off-axis 

contribution 

Total 

Exact 

Crude  S-P 

25 

8. 87698 (-1) 

8.93627 (-3) 

9. 77061 (-2) 

9.75652(-2) 

9. 260 (-2) 

50 

0.228548 

1.32791(-3) 

0.229876 

0.229878 

0.2324 

100 

0.560734 

-3. 80191 (-4) 

0.560354 

0.560418 

0.6189 

150 

0.801492 

2. 28596 (-5) 

0.801515 

0.801587 

0.8156 

200 

1-7. 64312 (-2) 

-1.88311  (-7)  1- 

-7.64314 (-2) 

1-7. 63845 (-2) 

1-7. 375 (-2) 

Np  -  18, 

G  »  72 

250 

1.50326 (-2) 

1.09488(-5) 

1.50435 (-2) 

1.50203 (-2) 

1.539(-2) 

500 

8. 17995 (-2) 

-4. 29986 (-4) 

8. 13695 (-2) 

8.13251(-2) 

8. 292 (-2) 

750 

0.208630 

4.65397(-4) 

0.209096 

0.209060 

0.2101 

1000 

0.369655 

-2. 76864 (-4) 

0.369378 

0.369373 

0.3692 

1500 

0.672094 

-3.98353(-5) 

0.672054 

0.672103 

0.7066 

2000 

0.860206 

-2.67366(-6) 

0.860204 

0.860256 

0.8678 

FIGURE  CAPTIONS 


Fig.  1.  Probability  distributions  (p^)  of  secondary  electrons  when  the  primary 
electrons  have  negative-binomial  distributions  (solid  lines)  or  arise  from  in¬ 
cident  light  with  a  Lorentz  spectral  density  (dashed  lines);  *  2,  G  »  6. 

The  curves  are  Indexed  by  the  number  M  of  degrees  of  freedom.  The  curve  marked 
*«>'  represents  the  Neyman  Type-A  distribution. 

Fig.  2.  Probability  distributions  {p^}  of  secondary  electrons  when  the  primary 
electrons  have  negative-binomial  distributions  (solid  lines)  or  arise  from  in¬ 
cident  light  with  a  Lorentz  spectral  density  (dashed  lines);  «  6,  G  2. 

The  curves  are  Indexed  by  the  nximber  M  of  degrees  of  freedom.  The  curve  marked 
represents  the  Neyman  Type-A  distribution. 

Fig.  3.  Paths  of  steepest  descent  of  the  integrands  in  (1.7)  and  (1.8)  when  the 
cumulative  Neyman  Type-A  distribution  is  being  calculated  for  N^  «  G  10.  The 
left-hand  curve  is  for  n  ■  75,  the  right-hand  for  n  ■  150.  The  small  circles 
indicate  the  saddlepolnts  of  the  Integrand. 

Fig.  4.  Cumulative  distributions  of  the  number  of  secondary  electrons  when  the 
primary  electrons  have  a  negative-binomial  distribution;  N^  5,  G  «  20.  The 
curves  are  Indexed  by  the  number  M  of  degrees  of  freedom,  the  curve  marked  *<x>' 
representing  the  Neyman  Type-A  distribution.  The  small  circles  denote  values 
of  the  cumulative  distribution  when  the  incident  light  has  a  Lorentz  spectral 
density  with  the  same  number  of  degrees  of  freedom. 

Fig.  5.  Cumulative  distributions  of  the  number  of  secondary  electrons  when  the 
primary  electrons  have  a  negative-binomial  distribution;  N^  ■  20,  G  >  5.  The  . 
curves  are  indexed  by  the  number  M  of  degrees  of  freedom,  the  curve  marked 
representing  the  Neyman  Type-A  distribution.  The  small  circles  denote  values 


of  the  cumulative  distribution  when  the  incident  light  has  a  Lorentz  spectral 
density  with  the  same  number  of  degrees  of  freedom. 

Fig.  6.  Cumulative  distributions  of  the  number  of  secondary  electrons  when  the 

primary  electrons  have  a  Ttegatlve-blnomlal  distribution;  N  -  18,  G  =  72.  The 

P 

curves  are  Indexed  by  the  number  M  of  degrees  of  freedom,  the  curve  marked 
representing  the  Neyman  Type-A  distribution.  The  small  circles  denote  values 
of  the  cumulative  distribution  when  the  incident  light  has  a  Lorentz  spectral 
density  with  the  same  number  of  degrees  of  freedom. 

Fig.  7.  Cumulative  distributions  of  the  number  of  secondary  electrons  when  the 
primary  electrons  have  a  negative-binomial  distribution;  •=  72,  G  =  18.  The 
curves  are  Indexed  by  the  number  M  of  degrees  of  freedom,  the  curve  marked  '<»' 
representing  the  Neyman  Type-A  distribution.  The  small  circles  denote  values 
of  the  cumulative  distribution  when  the  incident  ll^t  has  a  Lorentz  spectral 
density  with  the  same  number  of  degrees  of  freedom. 
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